Theory of Plasmon-assisted Transmission of Entangled Photons 
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The recent surface plasmon entanglement experiment [E. Altewischer et al., Nature (London) 
418, 304 (2002)] is theoretically analyzed. The entanglement preservation upon transmission in 
the non-focused case is found to provide information about the interaction of the biphoton and the 
metallic film. The entanglement degradation in the focused case is explained in the framework of 
a fully multimode model. This phenomenon is a consequence of the polarization-selective filtering 
behavior of the metallic nanostructured film. 
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Entanglement is one of the most strange prop- 

erties of quantum mechanics. Despite its puzzling charac- 
ter, this property, which is directly linked to the non-local 
nature of the theory, has been tested many times. The 
first experiments were performed with simple systems 
comprising two photons 0,01 or two atoms 0]. Recently, 
entanglement has started to be considered as a resource 
for diverse applications in quantum information theory. 
This has driven the interest in the demonstration of en- 
tanglement for systems involving many particles. Large 
systems are more prone to decoherence processes and, 
therefore, entanglement should be a very fragile property 
for them. For this reason the experiment of Julsgaard et 
al. 13] showing entanglement between two gaseous cae- 
sium samples, and the plasmon-assisted transmission of 
entangled photons shown by Altewischer et al. 0| have 
both attracted quite a lot of interest. This Letter is de- 
voted to Altewischer 's experience. Although some the- 
oretical aspects of the experiment have been treated in 
the original paper and in Ref. , a general (multimode) 
theory is still lacking. Here a complete detailed theory is 
derived and a thorough analysis addressing all aspects of 
the experiment is presented. 

In the experiment, pairs of correlated photons are gen- 
erated by spontaneous parametric down conversion. The 
desired input polarization-entangled biphoton state is 
obtained by appropriate manipulation of the generated 
photon pairs [lOj. This state is a quasi- monochromatic 
(A ~ 813 nm) quasi-plane wave (propagating along the Z 
axis) in the polarization singlet (1X1^2) ~ \Y1X2)) / 
where X and Y denote horizontal and vertical polariza- 
tion, respectively, and the subscripts {i = 1,2) label the 
first and second photon. After travelling along their re- 
spective trajectories parallel to the Z axis, each photon 
traverses a polarizer and is measured by a detector 
Dj [Pi is orthogonal to the Z axis and rotatable around 
it (the angle between the optic axis of P,; and the X 



axis will be denoted as Pi)]. These measured signals are 
electronically combined to obtain the rate of coincident 
photon detections. Such a setup allows one to determine 
the biphoton fringe visibility V^^ for fixed /32, which is a 
measure of the photons entanglement degree [l(] | . 

Photon 2 propagates freely from the source to P2, 
whereas photon 1 traverses a 200 nm thick Au film de- 
posited on top of a 0.5 mm glass substrate (Fig.^. The 
metallic film is drilled with cylindrical holes (200 nm di- 
ameter) arranged as a square lattice (period 700 nm). 
The transmission of photon 1 through the metallic film 
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FIG. 1: Optical elements in the trajectory of photon 1. The 
lenses, that constitute a confocal telescope, are only present 
in some parts of the experiment. The incident photon is 
represented by a monochromatic plane wave with wave vec- 
tor ki = {^.l,klz) ~ (0,fc). The first lens produces a fo- 
cused beam including many plane waves with wave vectors 
k2 ~ (q2,fc2z). The second lens gives rise again to a bundle 
of modes with ks = (q3,A;3z). The unit vectors P3, S3 indi- 
cate the p and s polarization directions, respectively, whereas 

P3 paraxi 
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(at A ~ 813 nm) is due to the phenomenon of extraor- 
dinary light transmission mediated by surface plasmon 
modes The chosen wavelength corresponds (al- 

most) to a transmission resonance ascribed to a surface 
mode at the metal-glass interface and propagating along 
the diagonals of the hole array, i.e., a (±1, ±1) mode. In 
order to investigate the effect of focusing on the entangle- 
ment behavior upon transmission, the film is positioned 
at the focus of a confocal telescope in some parts of the 
experiment (lenses' focal length / = 15 mm, telescope's 
numerical aperture 0.13). Altewischer et al. report en- 
tanglement preservation upon transmission when no tele- 
scope is used. The entanglement is however degraded 
when photon 1 is focused on the hole array and, more- 
over, the measured visibilities Vqo and V450 are different. 

A careful interaction model of the biphoton with the 
hole array should explicitly include in the wave function 
the quantum state of the solid. In a simplified monomode 
case (i.e., without telescope), the initial wave function is 
|$i„) = {\X1Y2) - \YiX2))/V2^ \S), where \S) is the ini- 
tial state of the solid. When the interaction has finished, 
the wave function j^out) can be written as: 

tx^xAXlY2) ® \S,x) +tY^xAYlY2) \Sy:,) 
- tx,YAXlX2)(g>\S^y)~tY,Y^\YlX2)(^\Syy), (1) 

where Iab are the transmission amplitudes for the var- 
ious channels, and \Sab) are normalized wave functions 
for the final state of the solid. The system's final quan- 
tum state is defined by post-selection, and for this reason 
|$out) only includes terms with exactly two photons. In 
other words, processes where, for instance, photon 1 is 
absorbed do indeed exist {toxi IOY2) ^ \Sox) — toYi 10X2) ® 
\Soy)), but they are not relevant for the visibility mea- 
surement because only coincident photons are registered. 
Notice that the final states of the solid must be taken into 
account, as may be clearly seen by considering the two 
following extreme situations: (?) all \Sab) are orthogonal 
to each other. In this case the solid and the biphoton 
can be entangled to a larger or lesser extent (depending 
on the tAB values) but the biphoton state (obtained by 
tracing over the solid) is always a mixture of factorizable 
states and it is therefore completely disentangled [l3| . 
Let us point out that in this case, after passage of photon 
1 , the solid incorporates a "which-polarization" informa- 
tion linked to the photons polarization state. This trans- 
lates into an entanglement loss, (ii) all \Sab) are equal. 
In this case the solid and the biphoton are completely 
disentangled, and the biphoton state can range from fac- 
torizable to maximally entangled depending on the Iab 
values. Since entanglement is preserved in some parts 
of Altewischer's experiment, {ii) is our model for the in- 
teraction process, i.e., the interaction does not introduce 
"which-way" labels in the solid. Such a theoretical frame- 
work explains why entanglement is preserved when the 
photon 1 is not focused. It is a simple consequence of 



two facts: first, no "which-way" labels are introduced in 
the solid, and second, the transfer matrix tAB for an or- 
thogonally incident plane wave (non-focused) on a square 
hole array is (by symmetry) proportional to the identity. 

Within the present model the visibility computation 
only requires the determination of the transfer matrix 
tAB for the employed optical set up. When the telescope 
is used (to focus the field at a spot on the film) , it is nec- 
essary to consider a multimode theory. The calculation 
proceeds as follows: (i) the electromagnetic field of pho- 
ton 1 is expanded in plane waves |q a) before and after 
each optical element (a = {p, s} denotes the two possi- 
ble polarizations; see Fig. {ii) Lens and film are de- 
scribed by their transfer matrices in the aforementioned 
basis, L(q2,qi), F(q2), respectively. {Hi) These matri- 
ces are combined to obtain the transfer matrix T(q3, q^^) 
of the telescope with the hole array inside it. {iv) The 
biphoton transfer matrix M[(q3,q2) for the whole set 
up (including polarizers) is then given by the tensor 
product of the transfer matrices for each photon, i.e., 
M(q3,qi) = Pi(q3)T(q3,qi) ® P2(0), where P, are the 
polarizers' transfer matrices. In the case of normal in- 
cidence, the telescope plus hole array transfer matrix 
T(q3,0) turns out to be: 

nh) J dq2e'^^^'»^"T^''^^'K-i(02)F(q2)M(02), 

(2) 

where R is the two-dimensional rotation matrix, n and 
A are the substrate refractive index and thickness, re- 
spectively (n = 1.52), k is the wave number, / is the 
focal length, and the remaining variables are explained 
in Fig. Note that the rapidly oscillating phase inside 
the integral means that the amplitude of |q3) is mainly 
given by the |q2) modes inside the telescope around the 
stationary phase condition q2 — nf(\^/{n ~ 1)A. 

To obtain Eq. Q a few approximations have been 
done. First, due to the low numerical aperture of the tele- 
scope, paraxial equations can be used. The lens transfer 
matrix L(q2,q]^) is given by: 

_^e^5t('J^-ii)'M(02)M-i(0i), (3) 

where the subscripts (j = 1, 2) refer to modes before and 
after the lens, respectively. Free propagation of a \<\2 (^2) 
mode inside the telescope along a distance z amounts to 
an extra phase that, apart from global factors, is given 
by exp(— zzq2/2fc) in the paraxial approximation. Sec- 
ond, the detectors are located at the back focal plane of 
an auxiliary lens placed after the polarizers. This im- 
plies that each point in the detector essentially collects 
one |q3) (in the limit of very large auxiliary lens aper- 
ture). For this reason the relative phases of the |q3) 
modes after the telescope are not relevant for the vis- 
ibility computation. Third, concerning the hole array, 
again due to the low numerical aperture of the telescope. 
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it is enough to keep the 0-th diffracted order (higher or- 
ders are not collected by the second lens), and therefore 
the film transfer matrix F(q) is q-diagonal. This matrix 
is numerically computed as explained in [l^ . Figure |21 
shows the computed transmittance when the hole array 
is illuminated by an orthogonally incident (non-focused) 
plane wave, and the film is tilted around the diagonal 
(compare to Figs, lb, Ic in 0). Despite the simulations 
do not exactly give the experimental peaks' heights and 
widths values, all main features are reproduced, including 
number and position of peaks, and overall order of mag- 
nitude. The behavior of the photonic bands as a function 
of parallel momentum is also correctly described. 
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FIG. 2: Hole array transmittance as a function of wavelength. 
The film is illuminated by an orthogonally incident plane wave 
(non- focused), and it is tilted (from 0° to 6°) around the di- 
agonal. The incident field linear polarization is perpendicular 
to the lattice diagonal in (a) and parallel to it in (b). For 
clarity the curves are offset 0.25%. 

The transfer matrix T(q3,0) corresponding to tele- 
scope plus film can be worked out analytically when 
q3 = 0, i.e., when the detector's aperture is extremely 
small. In this case the matrix before the integral in 
Eq. ^ disappears and the phase inside the integrand 
docs not depend on (j)2. If one writes down explicitly 
T(0, 0) and takes into account the symmetry proper- 
ties of the hole array, it can be shown that T(0, 0) is 
proportional to the identity [note that this result does 
not depend on the particular model for the numerical 
computation of F(q)]. For this reason the whole set up 
should again preserve entanglement when only the fol- 
lowing channels {q^ = 0) ^ (all q2) ^ (qa = 0) are 
considered. This means that the focusing of photon 1 
does not by itself degrade entanglement in every situa- 
tion but, rather, only when the transfer matrix Iab is not 
proportional to the identity. This could be easily checked 
in an experiment by inserting an iris before the detector. 

Figure 13 shows the visibilities obtained when the tele- 
scope is again employed, but now all channels [{q^ = 
0) (all q2) (all qg)] are considered. The visibility 
decreases as the telescope semiaperture increases (for 0° 



semiaperture the monomode case is obviously recovered) . 
The same behavior observed in |^ for the (±1, ±1) mode 
is reproduced by the simulations for A — 797 nm and 
telescope semiaperture 8°: V450 = 89%, Vo° = 37% (87 
and 73 in the experiment, respectively). The numerical 
discrepancy for Vqo can be essentially attributed to the 
fact that the computed transmission resonances are nar- 
rower than the measured ones, and the visibility is ver 
sensitive to this variable. It is also to be noted that in 
the employed wavelength is a bit larger than the resonant 
wavelength whereas we have computed the visibility at 
the resonance maximum itself (A = 797 nm). Our simula- 
tions show that Vqo grows for wavelengths larger than the 
resonant one (Vqo = 52% for A — 802 nm), whereas V450 
remains approximately the same. Note that the visibili- 
ties are not monotone functions of the semiaperture. This 
is due to the fact that for larger semiapertures, higher val- 
ues of 02 are included in the integral of Eq. and this 
permits the excitation of surface modes different from 
(±1,±1) (as can be indirectly seen in Fig.[21). 
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FIG. 3: Biphoton fringe visibility as a function of the tele- 
scope semiaperture. The visibilities are shown for two ori- 
entations of the second polarizer: /32 =0°, 45°. The chosen 
wavelengths are the resonances shown in Fig. |21 for no tilt: 
797 nm (corresponding to a surface mode propagating along 
the array diagonals), and 728 nm (corresponding to a mode 
propagating along the X or Y axes). 

The visibility reduction (as compared to the non- 
focused case) can be understood because the transfer ma- 
trix T(q3, 0) of the telescope plus film is not proportional 
to the identity anymore. This means that the system 
acts as a polarization-selective filter. The initial balance 
between the IX1I2) and II1X2) components (which is re- 
sponsible of the maximal entanglement of the input state) 
is therefore destroyed and as a consequence the entangle- 
ment is degraded. In the following it will be explained 
why do V450 and Vqo behave differently. 

Let us remind that, when Vp^ is measured, polarizer 
P2 is set with this P2 ■ To understand the behavior of Vg^ 
one can plot the field after the telescope due to a single 
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FIG. 4: Intensities (a, b) and polarizations (c, d) of mode Iqg) after the telescope as a function of {q3x,<l3y)- The excitation is 
one single linearly polarized photon (A = 797 nm). The incident photon polarizations are —45° in (a, c) and 90° in (b, d). The 
telescope semiaperture is 8°. The left color scale (arb. units) is shared by (a) and (b), and the right color scale corresponds 
to (c) and (d). In (c, d) the color scale represents the ratio of the minor to the major polarization elipse axes (for instance, 
±1 correspond to left and right circular polarization), and the small line segments indicate the orientation of the elipse major 
axis. The maxima (and minima) of all abscissa and ordinate scales correspond to a semiaperture of 63 — 0.1°. 



photon 1 incident with /32+90° linear polarization. This 
is so because the singlet biphoton state can be written as 
(1/32 /32 + 90°) - 1/32 + 90° p2))/V2 for every p2 (i.e., it is 
polarization isotropic). Vg^ will be 100% if there exists 
an orientation Pi for polarizer Pi that completely blocks 
the field due to photon 1 incident with linear polarization 
/32 + 90°. Otherwise, as /3i is varied, the collected inten- 
sity oscillates between a non-zero minimum and a maxi- 
mum, and Vfj^ < 100%. Such a polarization information 
is displayed in Fig. 0] Let us start with the analysis for 
V450. For —45° incident polarization of photon 1, only 
the (-1-1,-1) and (—1,-1-1) modes are excited. The out- 
put field is predominantly linearly polarized along —45° 
[Fig- Etc)], yielding high V450 (a low photon coincidence 
rate will be registered for /3i = -|-45°). On the other hand 
the analysis for Vb° goes as follows. For 90° photon 1 in- 
cident polarization, all (±1,±1) modes are excited, the 
output field is generally eliptically polarized [Fig. E^d)], 
and it includes various polarization directions. This fact 
is responsible for the decrease in Vqo (no /3i exists that 
nearly blocks the field). Note that the central region 
of both diagrams is linearly polarized along the incident 
polarization direction. This fits with the transfer matrix 
being approximately proportional to the identity (for this 
region) and with the 100% visibility expected for small 
detector aperture, as discussed previously. It is also to 
be noticed that Fig. Ejja) and (b) compare well with 
Figs. 2b, and 2d in |l4| [the stationary phase condi- 
tion mentioned above maps the diagram's side lengths 
shown here (26*3 = 0.2°) to the 26*2 = 17.6° value shown 
in Q; as opposed to circular fringes do not appear 
in Fig.0](a) and (b) in the range shown because multiple 
interferences in the substrate were not included in the 
calculation] . The A — 728 nm surface mode propagates 
along the (±1,0) and (0,±1) directions and, as a conse- 
quence of previous discussion, the roles of V450 and Vqo 
should be exchanged (this is indeed observed in the po- 
larization diagrams for A — 728 nm, not shown here for 



brevity). In fact, compared to A = 797 nm, the high and 
low visibility values are now exchanged, as it is distinctly 
seen in Fig. 13 where V450 is low and Vqo is high. 

In conclusion, a detailed multimode theoretical analy- 
sis of Altewischer et al. experiment has been presented. 
Entanglement preservation in the monomode case implies 
a particular model for the hole array-biphoton interac- 
tion, namely, this interaction cannot introduce "which- 
way" labels in the metallic film. Our model also repro- 
duces the measured results in the focused case. The en- 
tanglement degradation is understood as a consequence 
of the polarization-selective filtering behavior of the hole 
array for non-orthogonal incidence. A polarization anal- 
ysis explains the different values of the 0° and 45° visi- 
bihties. 



* Electronic address: l esteban.moreno@uam.es1 

[1] E. Schrodinger, Naturwissenschaften 23, 807 (1935). 
[2] E. Schrodinger, Naturwissenschaften 23, 823 (1935). 
[3] E. Schrodinger, Naturwissenschaften 23, 844 (1935). 
[4] S. J. Freedman and J. F. Clauser, Phys. Rev. Lett. 28, 
938 (1972). 

[5] A. Aspect, J. Dalibard, and G. Roger, Phys. Rev. Lett. 

49, 1804 (1982). 
[6] E. Hagley et al, Phys. Rev. Lett. 79, 1 (1997). 
[7] B. Julsgaard, A. Kozhekin, and E. S. Polzik, Nature 

(London) 413, 400 (2001). 
[8] E. Altewischer, M. P. van Exter, and J. P. Woerdman, 

Nature (London) 418, 304 (2002). 
[9] J. L. Velsen, J. Tworzydlo, and C. W. J. Beenakker, 
arXiv:quant-ph/0211103 (2003). 
[10] P. G. Kwiat et al, Phys. Rev. Lett. 75, 4337 (1995). 
[11] T. W. Ebbesen et al, Nature (London) 391, 667 (1998). 
[12] L. Martin-Moreno et al., Phys. Rev. Lett. 86, 1114 
(2001). 

[13] S. Hill and W. K. Wootters, Phys. Rev. Lett. 78, 5022 
(1997). 

[14] E. Altewischer, M. P. van Exter, and J. P. Woerdman, 



arXiv:physics.optics/0208033 (2002) . 



